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We study a mesoscopic ring with an in-line quantum dot threaded by an Aharonov-Bohm flux. 
Zero-point fluctuations of the electromagnetic environment capacitively coupled to the ring, with ui s 
spectral density, can suppress tunneling through the dot, resulting in a quantum phase transition 
from an unpolarized to a polarized phase. We show that robust signatures of such a transition can 
be found in the response of the persistent current in the ring to the external flux as well as to the 
bias between the dot and the arm. Particular attention is paid to the experimentally relevant cases 
of ohmic (s = 1) and subohmic (s = 1/2) noise. 
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The persistent current in a mesoscopic ring, penetrated 
by an Aharonov-Bohm flux, has been studied intensively 
in the past twenty years [lj. This equilibrium current is 
an indication of quantum coherent motion of charge car- 
riers in the ring. Recent studies show that the zero-point 
electromagnetic fluctuations in the leads can suppress 
quantum coherence, and effectively decrease the magni- 
tude of the persistent current . 

Theoretically, it is known that environmental dissipa- 
tion can even cause qualitative changes in the ground- 
state properties, i.e., it can drive the system across 
a quantum phase transition (QPT) into a dissipation- 
dominated phase. A popular model system is the spin- 
boson model 0, Q , describing a two- level "impurity" lin- 
early coupled to a bath of harmonic oscillators. Here, a 
QPT between a delocalized and a localized phase is found 
both for ohmic and subohmic damping |fj]. The more 
complicated Bose-Fermi Kondo model, with an addi- 
tional fermionic bath, also shows QPT between different 
ground states |E S @ ■ Several experiments have been 
suggested to observe such environment-induced QPT, 
utilizing sin gle- electron transistors with electromagnetic 
noise a, |lfj, , quantum dots coupled to ferromagnetic 
leads 3, or cold atoms in optical lattices [l2|. 

In this paper, we propose to use the persistent cur- 
rent in a metallic ring as a detector of a noise-induced 
QPT. A suitable setup Q consists of a ring with a small 
in-line quantum dot, capacitively coupled to an exter- 
nal circuit with dissipative impedance, Fig. 1. The key 
idea is that tunneling through the barriers can be effec- 
tively suppressed by charge fluctuations coupled to the 
ring electrons. Whereas Ref. studied the magnitude 
and the fluctuations of the persistent current in the pres- 
ence of ohmic noise, we focus here on the QPT caused 
by the dissipative environment. Modeling the external 
circuit by a RLC transmission line, we put emphasis on 
the cases of ohmic dissipation and a subohmic dissipa- 
tion with y/cJ spectrum density. They can be realized 
as LC-dominant and i?-dominant limit of the transmis- 
sion line, respectively I* 1 both cases we find robust 



signatures of a QPT between an unpolarized phase at 
small dissipation and a polarized phase at large dissipa- 
tion - in the latter, the charge state on the dot is dou- 
bly degenerate. Unlike previous proposals for detecting 
noise-induced QPT, our setup can be realized with cither 
ohmic or subohmic dissipation, and allows to assess the 
influence of finite temperature and bias. 

Model. The setup, shown in Fig. 1, consists of a 
polarizable environment coupled to the ring-dot sys- 
tem through the capacitors Cd and C a . For simplic- 
ity we assume spin-polarized electrons; adding the elec- 
tronic spin degree of freedom would enhance the mag- 
nitude of persistent current but not qualitatively 
change the physics studied here. Let us first consider 
the Hamiltonian for charges on the ring and the dot. In 
a small dot the level spacing, Ad, will be larger than 
the tunneling amplitudes t^ and t^. Regarding the level 
spacing in the arm, A a , we will separately analyze the 
cases of a small ring, A a 3> tL,tn, and a large ring, 
A a <C ii,i_R- For a small ring, low-temperature charge 
transport occurs only between the two topmost energy 
levels of the dot and the arm This leads to a two- 
level system i?tun = §<7z — Y ax > w ith e being the chem- 
ical potential difference between the dot and the arm, 
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FIG. 1: Ring with an in-line dot subject to a flux <E> and 
capacitively coupled to an external impedance Z. The dot is 
coupled to the left and right part of the arm through barriers 
with tunneling strength tz, and tR, and effective capacitors 
Cl and Cr. The impedance Z is modeled by an infinite RLC 
transmission line. 
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are Pauli matrices. 



The tunneling strength is 
2A£f,t a cos ip, where ip — 27r$/$ is 
the external flux in unit of the flux quantum $o = hc/e, 
and A = 1 (—1) for an odd (even) number of electrons 
on the ring. o~ z is related to the charge on the dot by 
Qd = v z e/2 + (N + 1/2) e, assuming N permanently oc- 
cupied states. In the opposite case of a large ring, the en- 
ergy spectrum in the arm is continuous, and progress can 
be made via bosonization of the arm electrons. This leads 
to a coupling between the pseudospin o~ z and the charge 
density fluctuations in the arm, i? h m = vp J2k>o kb\.bk — 
vfCz J2k>o \/ 7T k/2L(bk + b\) 0. It describes an ohmic 
bath J c h(w) = 27ra c hw6(w c — w) with a c h = 1/2. The 
cutoff is Lj c = VFk c , and vf the Fermi velocity. In the 
tunneling part, -f/tun, a factor y/ puj c /2iT now appears in 
the expression of A, where p is the density of states at 
Fermi energy in the arm. Thus, in both the small-ring 
and large- ring cases we arrive at a two- level system Ht un , 
but supplemented with an additional ohmic bath in the 
large-ring case. 

Now we consider the electromagnetic environment. 
The RLC transmission line, connecting to the dot and 
the arm, is described by the distributed resistances 
R, inductances L, and capacitances C per unit length 
(Fig.l). A Hamiltonian description of the environment 
can be obtained for an effective LC transmission line 
(R = 0). Allowing for an arbitrary distribution of L ni 
C n (n = 1,2,...), the Hamiltonian reads |2( H cm — 

2/(2C n ) + (^-0„_i) 2 /(2in) 
= iu>L\ + (iuCi + iloL2 + [iuiC-2 + 
is the charge on the capacitors 
C a and C d , and C^ 1 = C' 1 + + (C L + Cr)- 1 . The 
operators 4> n and Q n (n > 0) are conjugate operators sat- 
isfying [4>ri,Qm] = S nm . The system-environment cou- 
pling originates from the Coulomb interactions, i.e., the 
charging energy of the capacitors Cl and Cr, H coup — 
e/[2 (Cl + Cr)]Qo<t z . Hence the full Hamiltonian for the 
small- ring case reads H = H tun + H CO up + flcm, and for 
the large-ring case H = H tun + H coup + H cm + H ohm . 

To study the ground-state properties of H , we use Wil- 
son's numerical renormalization group (NRG) method 
[lfij . NRG is a non-perturbative method for impurity 
problems and was recently extended to impurity prob- 
lems with bosonic bath We diagonalize the bath 
degrees of freedom in H em and transform H into the stan- 
dard form of the spin-boson model 



Q 2 /(2C ) + EZi 
The impedance is Z (uj) 
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with the spectrum J(lo) = ^J^i^f^i^ ~~ F° r the 
small-ring case, J(u>) is solely from the electromagnetic 
contribution J em (u>) evaluated below. For the large-ring 
case, J(u>) = J cm (w) + Jch(w) where J c h(w) is from 7J hm 
as discussed above. 



The electromagnetic contribution J em (w) is related to 
the free retarded Green's function of Qn as J cm (w) = 
~e 2 /(C L +C R ) 2 Im((Q \Q }}^. Since ((Q \Q )) U = 
— \/[iuiZ(w) + C ], we have 



J em (oj) = —aujReZt (u) + ir/) 



(2) 



Here Z t (ui) = l/[Z _1 (w) + iuCo] is the total impedance 
describing the effective capacitor Co in parallel with 
Z(w), and the coefficient a = [eQ /(Qz, + Cr)] 2 . For 
the discrete transmission line in Fig. 1, Z(oj) is a compli- 
cated function, which in the continuum limit simplifies 
to Z(w) = [(R + iuj^/iuC} 1 / 2 [13. In both cas es, the 
asymptotic behavior of J cm is J cm = a^J R/(2C)\/u> for 
u> < R/L, 1/Vw and J em = a^/L/Coj for 1/VZC > 
uj ^> R/L. The crossover energy scale R/L separates 
the low-energy subohmic -Juj part from the high-energy 
ohmic u) part. 

For large resistivity R (i?-dominant leads), the y/oJ 
part of the spectrum will dominate the dissipation. In 
the opposite case (LQ-dominant leads), R/L is small, 
and the linear to spectrum will be relevant in the ex- 
perimentally accessible range of energies and temper- 
atures. For either type of the leads, the role of the 
charge modes in the large-ring case, J c h(w) cx u>, is only 
to renormalize the dissipation arising from J em (w) and 
will not change the quantum critical behavior. We will 
therefore describe both small-ring and large-ring cases 
using a "pure" bath with a power-law spectral density 
J(u>) = 2n 'auj s Lo^r s Q(uj c — u>), where the exponent s is 
determined by J cm (u>). The effective dissipation strength 
a is from J em (w) in the small-ring case, and in the large- 
ring case will be further renormalized by J c h(w). tu c is 
the effective cutoff and taken as our energy unit. 

The circulating persistent current operator is I c = 
(CrI l - C l Ir)/(C l + Cr), where I L and Ir are op- 
erators for particle current through the left and right 
tunnel barrier, respectively 0|. Due to (a y ) = the 
average displacement current (Il+Ir) through the exter- 
nal loop vanishes, and the expectation value of I c reads 
/ - (I L ) = -(Ir) - I {a x ). I = -e/2dA/d<p is the 
persistent current in the dissipationless limit. 

Phase transition. For a given experimental configura- 
tion, the bath exponent s and the dissipation strength a 
are fixed. A can be tuned in the range — tR\, th +£r] 
by scanning <&. The bias e is tunable by a gate volt- 
age on the dot. For both the ohmic and subohmic spin- 
boson models, a quantum critical point A c (a) separates 
the localized (A < A c ) and delocalized (A > A c ) ground 
states I n terms of the charges on the dot, these are 
the "polarized" and "unpolarized" states, respectively. 
Note that (o~ z ) can be viewed as an order parameter for 
this transition, with (a z ) =/= in the polarized phase. 
Whereas in the subohmic case a A c exists for any a (with 
A c — > as a — > 0), in the ohmic case a > 1 is required 
for the existence of a localized phase. Experimentally, 




FIG. 2: (a) Rescaled persistent current as function of A, for 
S = 0.2, a = 0.03; S = 0.5, a = 0.15; S = 0.8, a = 0.5; 
and S = 1.0, a = 1.2 (from top to bottom). The arrows 
mark the positions of A c . (b) Rescaled persistent current as 
function of bias e, at 5 = 0.5, a — 0.15 for A = 0.3 (triangles), 
A = A c = 0.18463 (squares), and A = 0.05 (circles). 

one can cross the critical point by scanning provided 
that \t L - t R \ < A c < t L + t R . 

Results. Our results were obtained via the bosonic 
NRG for the spin-boson model (2), using NRG param- 
eters A = 2 (logarithmic discretization), M = 80 (kept 
states), and Nf, — 6 (local boson states) First we 

focus on the ground-state persistent current I = (7 C ) and 
its response to A. The A-response is related to dl($) / 9$ 
which can be directly observed in experiment 0] . With 
the periodic background Jo removed, I/Iq = (<x x ) in- 
creases monotonously as A increases and saturates in the 
large-A limit. (Note that (a x ) ^ in both phases of the 
spin-boson model.) For the subohmic bath, < s < 1, 
I/Iq oc c(cv)A in the small-A limit, as shown in Fig. 2(a). 
c(a) increases with decreasing a and diverges in the limit 
a — ► 0. For the ohmic case, s = 1, our data extrapolated 
to A = 1 are consistent with the perturbative results 
I/Iq oc A K (") where nia) = a/(l — a) for a < 1/2 and 
k(o) = 1 for a > 1/2 0,01 • 

Approaching the quantum critical point A = A c , a 
singular contribution arises: I/Iq = (I/Iq) c + c + /_|A — 
A c | e + (I//o)rcg(A — A c ). The exponent 9 can be related 
to the longitudinal susceptibility exponent 733 0, Eof . 
defined as d(a z }/de oc |A — A c |~ 733 , using hyperscaling 
(see below), with the result 8 = 733/s — 1. Since 9 > 1, 
except for s = 1/2, I/Iq near A c is always dominated by 
the linear term of regular contribution (7/7o) reg , masking 
the power law in the raw data. However, the susceptibil- 
ity X11 = da x /dA = d(I/Io)/dA shows clear singularity 
in the range 1/3 < s < 0.92 where 6 — 1 < 1. For the 
i?-dominant leads, s = 1/2, a kink appears in the I/Iq 
curve at A c , manifesting the exponent 9 — 1 [2(J, see 
Fig. 2(a). Correspondingly, there is a finite jump at A c 
in the xn(A) curve. 

The singular behavior at the QPT can be nicely ex- 
tracted from the persistent current as a function of the 
bias e. 1(e) is a decreasing function and 1(0) depends on 
both a and A. 1(e) has a distinct small-e behavior in the 
different phases. For A > A c , 1(e)/ 1(0) oc 1 - ce 2 with 
a coefficient c depending on A and a; in the decoupled 




FIG. 3: Susceptibility X13 as function of the rescaled tunnel- 
ing strength A/A c . (a) subohmic case: S = 0.2, a = 0.03 
(circles) and S = 0.5, a = 0.15 (squares). The solid lines are 
for e = and dashed lines for e = 0.01. (b) ohmic case: S = 1, 
a = 1.2. A c fa 0.0714. For bias e = 10 -4 (dots), e = 10~ 6 
(squares) and e = 10~ 8 (up triangles), respectively. 

limit c = 1/(2A 2 ). For A < A c , 1(e) has a finite non- 
universal slope at e = 0. At the critical point A = A C! a 
power law appears: l-7(e)/7(0) oc lei 1 /* 13 . Such distinct 
behaviors are shown in Fig. 2(b) for s = 1/2. For other 
s values, while the behavior for the delocalized and lo- 
calized phases is similar, the slope of the critical curve at 
e = may be zero, a finite value, or divergent, depending 
on whether £13 (s) is smaller than, equal to, or larger than 
unity. The transition point 613 = 1 is at s = 1/3 (Fig. 
4(a)). For s = 1/2, the exponent <5 13 = 3/2, leading to a 
diverging slope of the critical curve, Fig. 2(b). 

To characterize the response of the persistent cur- 
rent to the bias, we introduce a transverse susceptibil- 
ity X 13 = -d(a x )/de\ e=0 = -(l/l )dl(e)/de\ e=a . Not 
only can this quantity be measured from the linear re- 
sponse of the persistent current to the bias, it is also 
related to the flux induced capacitance Gj, = d(Qd)/d& 
through (7$ = — e/(27rc)/oXi3 E3- I n tne delocalized 
phase (a z ) = 0, meaning X13 = C$ = for A > A c . 
X13 is finite for A < A c , and obeys a power law when 
approaching the quantum critical point from this side, 
X13 oc (A c — A)~ 713 . 713 is the transverse susceptibility 
exponent. In Fig. 3(a), we show X13 as a function of 
A/A c for s = 0.2 and s = 0.5 (symbols with solid lines). 
Since 713 changes sign at s = 1/3, near the critical point 
X13 approaches zero for s < 1/3 while it diverges for 
s > 1/3. This explains the different behavior of s = 0.2 
and s = 0.5 curves in Fig. 3(a). 

We examined the above scenario under a finite bias and 
finite temperatures. As shown in Fig. 3(a) (symbols with 
dashed lines), a small bias e = 0.01 smears the power-law 
singularity of X13 and leads to a rounded peak structure 
in the xi3(A) curve. The deviation from the zero-bias 
curve is most prominent near A c . Turning on a finite 
temperature, X13 are further suppressed on both sides of 
A c . However, in the regime T < e, the peak structure of 
X13 curve is always a pronounced feature. This provides 
a signature of the quantum critical point which is robust 
against finite bias and finite temperature. 

For the LC-dominant leads, s = 1, a Kosterlitz- 
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FIG. 4: Critical exponents (a) Sis and (b) 713 (crosses) as 
function of s, obtained from NRG. Also shown are the hyper- 
scaling results (a) 733(1 + s)/[2(7 33 — s)] and (b) 1 — 733(1 — 
s) / (2s) (circles with dashed line) . The solid lines are the val- 
ues obtained from hyperscaling relation by setting 733 = 1 in 
the s < 1/2 regime |2fl|. 

Thoulcss (K-T) transition occurs at a finite A c only for 
a > 1. In this regime, 1 — 1(e)/ 1(0) oc ce for both A < A c 
and A = A c , consistent with Si3(s — 1) — 1 17]. The 
transverse susceptibility follows Xi3 <x 5(A — A c ), in- 
stead of the (A c — A) -1 divergence naively expected from 
713 — ► 1 as s — ► 1. In Fig. 3(b), we plot the xi3(A) 
curve for s = 1, a = 1.2 under different bias. The peak 
structure for finite bias evolves towards a <5-peak in the 
small-e limit. A finite bias shifts the peak position from 
A c towards larger A, and broadens the peak dramati- 
cally. These Xl3 curves under finite bias thus present 
an well-defined precursor of the true QPT at zero bias. 
For finite temperature, T < e, similar peaks are observed 
with suppressed peak height. 

Finally, in Fig. 4 we show the two critical exponents 
£13 and 713. A scaling ansatz for the free energy |5| al- 
lows to derive hyperscaling relations between critical ex- 
ponents, with the results ^13 = 733 (1 + s)/ [2(733 — s)] and 
713 = 1 — 733(1 — s)/(2s), also shown in Fig. 4 with 733 
calculated by NRG. The numerical consistency confirms 
the validity of hyperscaling in the subohmic spin-boson 
model in the entire regime of < s < 1 5]. We also 
compare the numerical data with the hyperscaling results 
obtained by setting 733 = 1 for s < 1/2. Within numer- 
ical errors, most pronounced at s = 1/2 due to logarith- 
mic corrections, our data s uppo rt the analytical result 
733 = 1 for < s < 1/2 HEt Note that at s = 1/3, 
5i3 passes one and 713 passes zero, leading to changes in 
the critical behavior of observables as discussed above. 

Let us estimate the parameters in a typical exper- 
iment. We take the surface plasma frequency tu p — 
10 14 s _1 as the cutoff energy and energy unit. The upper 
boundary of the tunable dimcnsionless A is of the order 
1CP 3 — 10 -2 [13) an d the lower boundary can be very 
small for symmetric ring-dot system. The dissipation 
strength a s=1 ~ 1(T 2 and a s=1/2 ~ 10" 3 -10" 2 [13. For 
s = 1/2 (R-dominant leads) this parameter space covers a 
sufficient region of the phase diagram to make the obser- 
vation of QPT feasible. For s = 1 (LC-dominant leads) 
the a is too small, but the characteristic impedance Zq 



may be increased significantly by optimizing the circuit 
design and selecting suitable lead materials. Taking into 
account the a c h = 1/2 dissipation from charge fluctua- 
tions in the large- ring limit, it is still possible to realize 
an ohmic dissipation source with a > 1, and hence to 
observe the proposed features of the K-T transition. For 
a micron-sized ring, an in-plane magnetic field of 1 — 10 
T will effectively suppress the spin degrees of freedom to 
realize the spinless electrons employed here. 

Summary. We have identified robust signatures of a 
QPT, driven by zero-temperature equilibrium environ- 
mental noise, in a mesoscopic metal ring. Suitable ob- 
servables, showing quantum critical behavior, can be de- 
rived from the persistent current through the ring. We 
propose that this system can serve as a tool for study- 
ing the noise-induced QPT and probing the low-energy 
electromagnetic fluctuations in a mesoscopic circuit. 
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